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Chiral chemical potential does not cause the sign problem in the Monte Carlo simulation of lattice 
QCD. Using the chiral chemical potential, we study the chiral magnetic effect in two-flavor full 
QCD. We show that a strong external magnetic field induces an electric current in a chirally 
imbalanced QCD matter. The qualitative feature of the induced current is consistent with an 
analytical prediction. 
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1. Introduction 

Lattice QCD simulation is a powerful tool for nonperturbative analysis of QCD. However, 
it breaks down at a finite baryon chemical potential because of the famous sign problem. The 
fermion determinant becomes complex and the naive Monte Carlo sampling is invalid. For a small 
chemical potential, many methods has been proposed. It is difficult to apply such methods to 
a large chemical potential. In some special cases, we can exactly avoid the sign problem. The 
famous examples are isospin chemical potential and two-color QCD. These cases are well studied 
both in phenomenological studies and in lattice simulations. We here consider another possibility, 
that is, chiral chemical potential Compared with the isospin chemical potential and two-color 
QCD, the chiral chemical potential has been less studied so far. 

The chiral chemical potential /I5 is defined in the Dirac operator as 

D(Hs) = 7^D^+m + n 5 Y4Y 5 (1.1) 

The chiral chemical potential generates a finite chiral charge. The chiral charge is an imbalance 
between the left-handed and right-handed fermion numbers. The important property of the chiral 
chemical potential is that it does not cause the sign problem. 

Using the chiral chemical potential, we analyze the chiral magnetic effect. The chiral mag- 
netic effect is an electric current induced by a strong magnetic field in a heavy-ion collision [gj. A 
noncentral collision of two heavy ions produces a very strong magnetic field perpendicular to the 
reaction plane. The strong magnetic field fixes the spin and momentum directions of the quarks 
depending on their chiralities. The positive-helicity particles and the negative-helicity particles 
generates electric currents in the opposite direction. If the left-handed and right-handed chiralities 
are symmetric, these two contributions exactly cancel out, and the net current is zero. If the chiral- 
ities are imbalanced, the net current is finite. In QCD, such an imbalance is locally generated by 
the axial anomaly and the topological fluctuation of the background gauge field. The chiral mag- 
netic effect is a direct experimental evidence of the topological fluctuation or the event-by-event 
CP violation g [|. 

In this study, we use the chiral chemical potential, instead of the topological fluctuation. This 
is a different approach from other lattice simulations for the chiral magnetic effect [§, 0, §]. The 
chiral chemical potential generates the chirally imbalanced QCD matter as an equilibrium state. 
Strictly speaking, the chiral charge is not conserved quantity, and this prescription is a kind of ap- 
proximation of real QCD. However, the chiral chemical potential is quite convenient for theoretical 
studies. The chiral chemical potential was introduced in a phenomenological work, and the induced 
current of the chiral magnetic effect was derived as 

j = 2tfte1 B ( L2 ) 
using the Dirac equation coupled with the background magnetic field [fj]. 

2. Simulation setup 

We performed the two-flavor Hybrid Monte Carlo simulation with the chiral chemical poten- 
tial. We used the Wilson gauge action and the Wilson fermion action. The color number is N c = 3. 
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For sake of simplicity, we considered two fermion flavors with the same mass m and charge q. This 
approximation simplifies the numerical simulation, especially the hybrid Monte Carlo algorithm. 
The Wilson-Dirac operator with the chiral chemical potential is 



[D w (^)]x,y = (l-ri)U i (x)S x+h + {l + yt)uj(x-t)8 x _ 



<■)■ 



— K 



[1 - Y 4e a ^)U 4 (x)8 x+ly + (1 + YAe- a ^)Ul(x-A)8 x _ ly 



e ±aw5 _ cosh(a/X5)±7 5 sinh(<3/x 5 ) . 



(2.1) 
(2.2) 



This is the simplest choice of the chiral chemical potential in lattice QCD. The form of the chiral 
chemical potential is analogous to baryon chemical potential ||]. The Dirac operator satisfies the 
relation 



750(^5)75=^(^5) 



(2.3) 



From this relation, we can immediately show that the fermion determinant detD(/i 5 ) is always 
positive real and thus the sign problem does not occur in the case of even flavors. 

For the analysis of the chiral magnetic effect, we introduced an external magnetic field. On 
the lattice, the U(l) gauge field is introduced as the Abelian phase factor Un(x). For the external 
magnetic field, the SU(iV c ) link variable is replaced as 



Un(x) ->■ u^(x)U^(x) 



(2.4) 



only in the Dirac operator (2.1). The kinetic term of the U(l) gauge field is not introduced in the 
Lagrangian. The strength of the magnetic field is an external parameter and is not affected by 
dynamical effects. For applying a constant magnetic field B in the X3 -direction, the phase factor is 
set as 



u\(x) = exp(—iaqBN s X2) forjci =aN s 

u%{x) = exp(iaqBxi) 

Uu (x) = 1 for other components 



(2.5) 



with a 2 qB = (2n/Nj) x (integer) 



3. Chiral magnetic effect 

In this section, the lattice gauge coupling is jS = 2N c /g 2 = 5.32144 and the hopping parameter 
is K = 0.1665. These values correspond to the lattice spacing a ~ 0.13 fm and the pion mass 
m % ~ 0.4 GeV [[Hj]. 

In Fig. [IJ, we plot the chiral charge density 

n 5 = -a 3 (\j/Y475¥) = a 3 (YlYL - Wr¥«) ( 3j ) 

scaled by the lattice unit. At a finite chiral chemical potential, the chiral charge density is finite, 
i.e., the system is chirally imbalanced. The chiral charge density increases as the chiral chemical 
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Figure 1: The chiral charge density ns. The lattice sizes are Ng xN t = 12 3 x 4, 12 3 x 8, and 12 



potential increases. The qualitative behaviors are different between the confinement phase (N t = 8 
and 12) and the deconfinement phase (N t = 4). Note, however, that the chiral charge density 
saturates in afc > 1.0. This artificial behavior is called as saturation, which is also known in 



the cases of isospin chemical potential and two-color QCD flUQ. The lattice calculation works only 
below the saturation, i.e., a/15 < 1.0 in the present setup. 

For analyzing the chiral magnetic effect, we applied the external magnetic field to this chirally 
imbalanced QCD matter. We calculated the vector current density 



(3.2) 



We applied the magnetic field in the X3-direction and measured the transverse component j\ and 
the longitudinal component 73. The two transverse components are the same, j\ = }%, because of 
the rotational symmetry. The simulation was done in the deconfinement phase (N t = 4), where the 
chiral magnetic effect is expected in heavy-ion collisions. 

In Fig. ^| the transverse component j\ is plotted as a function of the magnetic field B (left) 
and of the chiral chemical potential /I5 (right). The transverse component is always zero because it 
is irrelevant for the chiral magnetic effect. As shown in Fig. ||, the longitudinal component j'3 is a 
linearly increasing function of B at finite /I5 and a linearly increasing function of at finite B. We 
can parametrize the induced current as 



j 3 = cfCNdatHsqB 



(3.3) 



The factor AW = N c X Nf = 6 is the number of quarks with the same charge. This functional form 



is consistent with the analytical approach ( |1.2| ). The overall constant C is 0.013 ±0.001 in the 
present lattice simulation. To compare this value with the analytical approach, we need to calculate 



the renormalization constant because the local vector current (3.2) is not renormalization-group 
invariant. When we carefully estimate the renormalization effect and systematic errors, we can 



evaluate QCD corrections to the analytical formula (1.2). A possible correction has been suggested 
in a phenomenological work [|T2|]. 
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Figure 2: The transverse component j\ of the vector current density. The left and right panels are plotted 
as a function of qB and /X5, respectively. The lattice size is Nj x N t = 12 3 x 4. 




Figure 3: The longitudinal component j$ of the vector current density. The notation is the same as in Fig. 2. 
The black dashed line is a linear function (3.3). 



In this study, we have succeeded in observing a finite induced current. This is completely 
different from other lattice simulations. In the standard lattice QCD without the chiral chemical 
potential, the induced current of the chiral magnetic effect cannot be observed [g]. Since lattice 
QCD can reproduce nontrivial topological sectors, one might think that it can reproduce the induced 
current. However, this is not so easy. In the usual lattice simulation, the number of the topological 
charge is typically 0(1), and the lattice volume is typically a 4 V ~ 0(1O 5 ). The topological charge 
per volume is thus 0(1O~ 5 ). This is rather small. In our simulation at a finite chiral chemical 
potential, the chiral charge density is independent of the volume and its value is O(10 _1 ). Owing 
to such a large value, we can directly observe the induced current of the chiral magnetic effect. 
This is a great advantage of the chiral chemical potential. 



4. Confinement/deconfinement phase transition 

We also analyzed the phase structure of a chirally imbalanced QCD matter. The phase structure 
in H5-T plane was studied in phenomenological models [13, 14, 15]. We examined the temperature 
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Figure 4: The expectation value of the Polyakov loop P at ajXs = 1.0. The lattice sizes are iV 3 x N, = 
8 3 x 4, 12 3 x 4, and 16 3 x 4. 

dependence of the Polyakov loop by varying the lattice gauge coupling /3 . The Polyakov loop in 
full QCD is a phenomenological criterion for confinement/deconfinement. The hopping parameter 
is fixed at K = 0.1665. The magnetic field is not applied here. 

In Fig. I], we plot the Polyakov loop at a/15 = 10. The Polyakov loop rapidly rises at /3 ~ 5.27. 
This corresponds to a confinement/deconfinement phase transition. The low-temperature side is the 
confinement phase, and the high-temperature side is the deconfinement phase. We determined the 
order of the phase transition from the dependence on the spatial volume V = a 3 N 3 . As shown in 
the figure, the Polyakov loop is almost independent of the spatial volume. We checked that its 
susceptibility is also independent of the volume. This scaling behavior suggests that this transition 
is a crossover. We also calculated at a}!*, = and 0.5, and found that the situation is the same. Thus, 
we conclude that the order of the confinement/deconfinement phase transition does not change in 
< aji.$ < 1.0 in the present setup. Since the nature of the phase transition depends on the quark 
mass, the calculation at a different quark mass can leads a qualitatively different result, such as a 
first-order phase transition or a critical endpoint. 
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